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1. Theoret i cal an d Numerical Analyses 

The problem of vacillation in a baroclinically unstable flow field is being studied 
as a task for this grant through the time evolution of a single nonlinearly unstable 
wave. To this end a computer code is being developed to solve numerically for the 
time evolution of the amplitude of such a wave. The final working code will be the 
end product resulting from the development of a hierarchy of codes with increasing 
complexity. The first code in this scries has already been completed and is presently 
undergoing several diagnostic analyses to verify its validity. The development of this 
first code is detailed below. 

The first model to be studied here is what is commonly known as the ” Eady 
problem". This problem is concerned with the baroclinic instablity of an infinite 
fluid on an f-plane and without, potential vort.icity. The basic state (low for this 
model is for constant shear. IJy assuming the potential vorticity to be initially zero 
it must then remain so for all time t 0. With these assumptions the equation 
governing the pressure perturbation, p, for the flow in the interior of the layer takes 
the following form: 
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It was decided to begin the hierarchy of codes by first tackling the Eady prob- 
lem propper, i.e. the baroclinic instability problem without Ekrnan dissipation. 
This is done with the knowledge that such a model might not develop the vacilla- 
tion phenomenon being investigated. However, by supressing the Ekrnan damping 
graeter simplicity in coding is achieved without altering any of the fundamentals of 
the numerical techniques employed in the computer code. Thus, with these consid- 


erations in mind Equation (l) was then solved subject to the following boundary 
conditions: 
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Note, that in this case the governing equation, (l), is linear while the upper 
and lower boundary conditions, (3), are nonlinear. This fact then allows for the 
pressure perturbation function to be expressed in the foilwing way: 
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where y. - (lc 2 ~ l 2 n 2 )S | l/2 , and the coefficients {Ak,hBk.i) are equal to (A' 

B'-k . (H° insure the pressure perturbation, p is real. Note that the vertical struc- 
ture of each component of expression (4) depends on the horizontal wavenumber. 
( k 2 — l 2 i r 2 ) l/ 2 . and on the stratification parameter. S. The amplitude coefficients. 
Ak,i and Bk.i . are functions of time and their evolution from some prescribed initial 
conditions is the subject of this study. 

Upon substituting the solution form (4) into the boundary conditions (3) and 
performing the necessary algebraic manipulations, an infinite set of nonlinear, or- 
dinary differential equations in time for the amplitude functions results. Note, that 
the solution form (4) satisfies the governing equation exactly. These equations take 
the following form: 
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where FV/ and Gk.i arc complicated nonlinear functions of all the amplitude co- 
efficients Ak,i and Bk,i . Clearly the resulting infinite systemof ordinary differen- 
tial equations must be truncated to some finite value of K, ( - K < k < K) and 
L,( — L < l < L) in order to obtain a numerical solution to the original system. 
The code developed for this period allows for the solution for any arbitrary value 
of both K and L and which need not be the same. 

The solution representation of the pressure perturbation function, (4), lends 
itself to the spectral form of numerical solutions to partial differential equations. 
Consequently this will be the numerical method employed in the development of 
the code. Furthermore, the pseudo-spectral method is used in the code developed 
for this problem. This means that the products and sums needed for generating 
the functions Fk.i and Gk.t were performed in the physical space, (x, y), while the 
result is obtained in the wavenumber space, (kj). A computer code for solving the 
coupled truncated set of ordinary differential equations was written. The integra- 
tion is started at t = 0 with some known intial values for all Ak.i and fik.l- the 
forward integration in time is accomplished via a second order Adams- Bashforth 
time stepping scheme. The basic algorithm followed in the code is the following. 
First, a value for both K and I. is chosen. Knowing the initial values for all the Ak.i 
and Bk.i\ p and its spatial derivatives are then computed from the expansion ( I) 
at the physical grid points. The grid points multiplications and additions are then 
calculated and Fourier transformed to obtain the right hand sides of Equations (5). 
Finally, the integration is stepped once in time and the new values of all Ak.i and 
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Bkj at the new time are obtained. The process is then repeated untill the final 
time is reached. 

This code is undergoing several diagnostic proceedures at the present time. All 
of the subroutines have been independently verified and checked. However, when the 
code was run in its entirety a stiffness problem was encountered. This problem was 
manifested by having overflow errors whenever the final integration time exceeded a 
specific value. This value was found to depend also on the number of terms retained 
in the expansion (-1), i.e. the value of both K and f,. The longest final time that we 
were able to run the code to was Id days and that was for values of K and L both 
equal to 5. Figures 1 and 2 show the time evolution of the amplitudes and B\ \ 
belonging to a single linearly unstable wave and for initial values of .015 and .023 
respectively. These figures show that after approximatly nine days the amplitudes 
of the higher harmonics for both components of the wave begin to exert substantial 
influence. However, these amplitudes do not show any vacillation tendencies for the 
time period shown. This is expected since it is well established that the period of 
vacillation brought about through barodinic instablity is roughly between two and 
three weeks. Thus in order to observe vacillation in our analysis integration to a 
much longer final time is necessary. Certainly the integration time of If days is 
nowhere sufficient. 

This problem of integration blow-up for larger times has not been resolved yet, 
although quite considerable effort have been expended towards resolving it. The 
plan for the second year of funding is to first resolve this problem satisfactorily and 
then to introduce F.kman damping conditions on the upper and lower boundaries. 
It should be kept in mind, however, that experience from annulus experiments 
suggest that vacillation occurs only when the upper lid of the annulus is removed. 
Thus if vacillation is not found with Fkman damping on both boundaries then that 
condition will be relaxed on the upper boundary alone. 
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2. Atomspheric Data An alysis 


The question of whether a vacillation phonomenon, similar to that observed 
in the annulus, exists in the atmosphere has not yet been satisfactorily answered. 
There is reasonable evidence on the existance of a two to three week cycle in the 
zonal index calculations at mid-latitudes in both the northern and southern hemi- 
spheres. However, the origin of these cycles, i.e. whether barotropicaliy or baro- 
clinically induced is still being investigated. An important task undertaken in the 
present contract is to address the existance and the origin of the 2-3 week cycle in 
the atmosphere. The originaaly proposed task called for the study of this oscillation 
phenomenon through analyzing available zonal indicies. 

During this first year of funding an agreement has been established between 
the P.I. and the Laboratory for Atmospheres, (GLA), at the Goddard Sapce Flight 
Center to use the FGGE data sets available at the GLA. With this agreement 
the GLA is providing for the use of all tapes containing the needed data as well 
as supplying the P.I. with the necessary software and compuetr time needed to 
complete the proposed analyses. As an initial step in gaining sulhcient working 
fluency in handling the FGGE data sets it was decided to first experiment with 
a small portion of these data. This portion is the first Special Observing Period, 
(SOP 1), which covers the period between Jan. 1 - Feb. 28. 1979. This data set 
is gridded and well organized making it convenient for the specific use required 
here. Consequently, the necessary computer codes required for generating the zonal 
indicies were written and verified. 

In generating the zonal index the established conventioal methods were fol- 
lowed. Briefly, the zonal index, R, in this work will be defined in the following 
manner: 
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where A denotes longitude, 4 > latitude and V the horizontal velocity vector at a 
given grid point and pressure surface with u and v representing the eastward and 
northward components respectively. 

Figures 3 and 4 show the zonal index variation at the 500mb surface over a 
period of 40 days starting with Jan. 5,1979 and calculated for latitude positions of 
26 degrees north and 38 degrees north respectively. The time perod shown in these 
figures is appropriate for the SOP l experiment. Note, that since the level Illb of 
the FGGE data are given at 6 hours intervals the records shown in these figures are 
plotted for that same interval in time. It is possible to detect a 2-3 week cycle in 
this time record as well as several higher frequency cycles. These figures are shown 
here specifically for demonstration purposes and no further comments will be made 
with regards to thier meaning or utility. 

For the next funding year zonal indicies for the entire FGGE year will be 
generated at several latitude positions covering the region between .30 and 60 degrees 
of latitude for both the northern and southern hemispheres. These indicies will then 
be subjected to several numerical filters to establish the occurance of the 2,-3 week 
period cycles. For the time periods exhibiting the strongest signal of such a cycle 
full energy budgets will be evaluated and analyzed to help in establishing the origin 
of such cycles and their relationship to the vacillation phnomenon in question. 
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FIGURE 4 


